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ABSTRACT 

Larval transport and oceanographic conditions experienced by pelagic larvae in 
the Northwestern Hawaiian Islands were simulated using an individual-based approach 
to track daily movements in a Lagrangian modeling framework. These advection- 
diffusion models were configured with l°xl° resolution, monthly geostrophic currents 
estimated from satellite altimetry. Larval dispersal was simulated for each month of the 
year from 1993-2002 for 3-, 6-, and 12-month laiwal durations. Four release locations 
were evaluated: Midway Island, Maro Reef, Necker Island, and Oahu. Larval retention 
was evaluated by tabulating successfully simulated settlement, which was scored based 
on larval proximity to release sites after completion of the pelagic duration. Sea surface 
temperature and chlorophyll concentration at each daily larval location were tabulated 
utilizing similar resolution, satellite remotely sensed data products (NOAA Pathfinder 
AVHRR SST and SeaWiFS ocean color), and these in situ values were integrated over the 
entire larval duration for each larval track. These oceanographic variables are of critical 
importance in the early life histoiy because of their hypothesized relationships to larval 
growth and feeding success, both critical determinants of laiwal suiwival and successful 
recruitment. The sea surface temperature and chlorophyll histories experienced by 
successfully settling larvae display strong seasonal and interannual patterns, which were 
decomposed using generalized additive models (GAMs). These patterns may be useful 
towards understanding episodic recruitment events, as well as for posing hypotheses 
towards understanding the mechanisms underlying spawning seasonality. These transport 
dynamics and oceanographic patterns have general implications for a variety of vertebrate 
and invertebrate metapopulations in the Northwestern Hawaiian Islands and their 
effective management. 
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INTRODUCTION 

Temporal patterns of reproduction are a widespread phenomenon in both plant 
and animal ecology. Organisms can maximize their fitness by propagating at times which 
are optimal for maximum reproductive output and/or enhanced suiwival of their young. 

It is commonly thought that the latter is of more importance for highly fecund aquatic 
species which broadcast their abundant young into the pelagic environment (Johannes, 
1978; Thresher, 1984). Several scales of temporal variability may be of importance in 
the timing of reproduction. Diel patterns, such as spawning near dawn or dusk, may 
be important to minimize predation on both the spawning individuals and their pelagic 
propagules (e.g., Doherty, 1983). Lunar patterns, such as spawning near spring tides 
(full and new moon) may be related to key variables which change consistently on a 
monthly scale such as tidal currents and moonlight illumination (e.g.. May et al., 1979). 
Seasonal patterns, such as spring or summer spawning, may be related to key variables 
which change consistently on an annual scale such as currents, plankton blooms, and 
temperature (Johannes, 1981). Seasonality in spawning has been well documented in a 
variety of Hawaiian fish and invertebrate species (e.g., Itano, 2000; Lobel, 1989; Randall, 
1961; Reese, 1968; Walsh, 1987). Various hypotheses have been put forth to explain 
such seasonality. Johannes (1978) has argued that predatory losses on pelagic propagules 
have been a driving selective force for spawning seasonality. More recent views of 
pelagic laiwal transport have emphasized the importance of retention with the view that 
many marine populations are more closed than open (Cowen, 2002; Cowen et al., 2003; 
Jones et al., 1999; Kingsford et al., 2002; Leis, 2002; Mora and Sale, 2002; Robertson, 
2001; Sponaugle et al., 2002; Swearer et al., 1999). While predation and retention issues 
may be important, the predominately oligotrophic pelagic environment has led some to 
suggest that larval food supply is the single most important factor governing the numbers 
of marine fish (Cushing, 1972). Reese (1968) suggested that the different spawning 
seasons used by ecologically similar species of hemiit crabs were a mechanism to 
reduce competition for pelagic larval food, and that there would be selective advantage 
to offsetting reproductive periods if larval food supply were a limiting factor. Larval 
food supply involves spatial and temporal patchiness, and the species composition 
of the phytoplankton and microzooplankton is critically Important (Lasker, 1975). In 
addition to starvation issues, variability in food supply has been shown to be an important 
determinant of larval growth and subsequent survival (e.g.. Booth and Alquezar, 2002). 
Faster growth has been hypothesized to favor survival by reducing cumulative predatory 
mortality (e.g., Anderson, 1988). Leis and Carson-Ewart (1999) suggest that larger 
size is an important factor for evading predation during the settlement process, citing 
fin erection behavior and importance of speed when fleeing predators, based on field 
experiments with coral trout larvae. It is possible that small size may be an advantage 
for remaining undetected during settlement; however, the advantages of being larger 
in the plankton probably outweigh the disadvantages, considering the gauntlet of size- 
based pelagic predators (reviewed by Purcell and Arai, 2001; Zaret, 1980). Additionally, 
larvae which grow faster may retain a size and survivorship advantage during the critical 
first few weeks post-settlement on the reef (Bergenius et al., 2002; Booth and Hixon 
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1999; Sponauglc and Pinkard, 2004). While size is most directly a function ofage, 
both food and ambient temperature have been shown to ha\ c a strong positive clVcct 
on larval growth (c.g., Buckley ct al., 2004). Clearly, in addition to physical retention, 
there is a suite of other considerations critical in the early life history sur\ ival of insular 
species. These important ecological considerations can be synthesized within a computer 
simulation using available tools and data products. 

Earlier works have used advection-diffusion models to examine lar\ al transport 
and retention (e.g., Griffin et al., 2001; Hill, 1991; Polovina ct al., 1999; Siegel ct al.. 
2003). Few such applications have integrated the oceanographic conditions experienced 
by individual larvae directly into the model. With the availability of remotely sensed 
data products, it is logical to incorporate these environmental fields into the computer 
simulation framework, particularly with individuabbased modeling approaches (c.g., 
Million et al., 2002). Sea surface temperature (SST) and chlorophylba concentration 
are widely available from a variety of satellite sensors, and both of these variables may 
have important linkages to the ecology of early life histor>^ stages, as described above for 
growth and mortality. The goals of this paper are to examine, via computer simulation 
and use of remotely sensed environmental data, the seasonal and interannual components 
of laiwal retention, transport, growth, and sur\4val in selected regions of the Hawaiian 
Archipelago. 



Simulated lar\^al releases were stratified by year (;?=10: 1993-2002), month (;7=12: 
January-December), locations spanning the Hawaiian Archipelago (/7=4: Midway, Maro, 
Necker, and Oahu, see Fig. 1), and larval duration (/?=3; 3, 6, and 12 months) to yield 
a total of 1,440 model treatments. These larval durations were chosen to bracket the 
known durations of several commercially important species of lobster and deep-water 
bottomfish. Five thousand simulated lar\^ae were released for each model treatment for 
a total of 7.2 million individuals. Each individual was tracked daily for the entire lar\ al 
duration in Lagrangian fashion using the following equations: 



where x represents longitude, v represents latitude, t represents time in days, ii represents 
the EastAVest component of the current speed, v represents the Nortli/South component 

of the current speed, cos( v^) adjusts distance by latitude to account for the spherical 
coordinate system, and D is the diffusivity coefficient (500 mVsec following Polovina 
et al., 1999). The currents utilized in this study were monthly 1° latitude/longitude 
resolution geostrophic flow fields calculated from satellite altimetry obtained from 
CNES/AVISO/SSALTO (CLS Space Oceanography Division, France). Integrated SST 
and chlorophyll-a histories encountered daily by individual lar\^ae were tabulated daily 
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using inteipolations from monthly 1° latitude/longitude resolution data grids. SST data 
were obtained from the MCSST (NOAA Pathfinder AVHRR satellites) product from 
NASA/JPL. Chlorophyll-a data was obtained from the Sea-viewing Wide Field-of-View 
Sensor (SeaWiFS) instrument on board the Seastar satellite, integration was performed 
by averaging the daily SST or chlorophyll-a interpolations over the entire pelagic 
duration. Additionally, only averages from the subset of larvae scored as successfully 
being retained by the source site following completion of the entire pelagic duration 
were tabulated to each treatment. A 140-km radius for scoring larval retention was used, 
similar to Polovina et al. (1999). This is an arbitrary threshold utilized as a compromise 
to achieve some larval settlement at this level of propagule sample size (/7=5,000). 

Since settlement success was evaluated in a relative manner (e.g., comparing between 
years, months, sites, pelagic durations), this exact dimension is not critical. For each 
of the 1,440 treatments, the following were tabulated: the number of larvae scored as 
retained, the number of laiwae scored as settling at Oahu, the number of larvae scored 
as not settling at any of the four sites, the average SST encountered by this subset of 
retained larvae, and the average chlorophyll-a encountered by this subset of retained 
laiwae. SST was available for the entire temporal duration of this analysis; however, 
chlorophyll-a data were available only from 1997 onwards. The advection-diffusion 
model was written in the open-source software XBASIC (http://www.xbasic.org) and run 
on an Intel P4 Windows XP system. Generalized Additive Models (GAMs) were used 
to identify relationships between a suite of response variables (retention, SST history, 
or chlorophyll-a history) and a suite of predictor variables (year, month, location, and 
pelagic duration). GAM is a relatively new analytical technique (Hastie and Tibshirani, 
1990) which is useful when the predictor variables have unknown a priori and possibly 
nonlinear effects upon the response variable. GAM analysis was canded out using the 
analytical software package S-Plus v. 6. 1 .2r2 on an Intel P4 workstation using Redhat 
LINUX 7.3 OS. Six GAM analyses were perfonned as outlined in Table 1, with each 
utilizing a different suite of predictor variables as described. The graphical output in the 
form of smoothing splines and comparative categorical effects serve as the primary basis 
for inteipretation, using standard error terms to indicate statistical significance. The GAM 
plots indicate the modeled relationships between the suite of predictor variables and the 
response variable, and the cumulative (hence additive) sum across each predictor function 
scaled by an intercept results in the predicted value for the response variable. Higher 
values along the y-axis indicate a higher contribution towards the predicted value for the 
response variable, and vice-versa. 



RESULTS 

Retention varied from a low of 0/5,000 (scored for 39 different treatments) to a 
high of 3,908/5,000 (scored for October, 1993, Midway, 3-month release). Retention was 
strongly related to year and larval duration, and had weak relationships to month and 
site (Figs. 2A-2D). Laiwal settlement at Oahu was strongly dependent on spawning site 
(Fig. 3). The Oahu settlement GAMs were run separately for the spatial effect and other 
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variables due to the numerous zeros in the data for sites farther from Oahu. For example, 
85% of the Midway runs yielded 0 larval settlement to Oahu, and the small amount of 
sueeessful settlement was aeeomplished only at the longest larval duration (Table 2). 

This laek of data eontrast in other sites effeetively weighted the GAM primarily towards 
the Oahu site, leading to diffieult interpretation. Henee, the GAM was run separately for 
eaeh site (Figs. 4-7), with the last GAM being a simple retention analysis for Oahu only 
Strong yearly effeets were observed in all sites, with weak monthly effeets, and duration 
only beeoming important at Neeker and Oahu. Oahu settlement was favored by a longer 
larval duration from Neeker (Fig. 6C, the elosest site to Oahu), and by a shorter larval 
duration from Oahu itself (Fig. 1C). 

Larval nonsettlement was eursorily examined in this analysis. Considering that 
there are abundant other sites available for larval settlement, this result should be treated 
with eaution. Flowever, by examining the larvae that did not settle at any of the four sites, 
some useful hypotheses ean be posed for further analyses. The data suggest that yearly 
and monthly effects may be relatively weak and that perhaps there is a spatial component 
involved with Oahu exhibiting higher rates of larval loss (Fig. 8C). As expected, a longer 
larval duration is positively correlated to larval non-settlement (Fig. 8D). 

SST and chlorophyll-a histories had strong yearly, monthly, and site relationships 
(Figs. 9-14), with the expected deterioration of a seasonal effect at a 12-month larval 
duration (Figs. 1 IB, 14B). For a 3-month larval duration, SST history was optimized by a 
July- August spawning, peaking in late July/early August (Fig. 9B). For a 6-month larval 
duration, the optimal spawning with respect to SST history is offset accordingly to May- 
July, peaking in June (Fig. lOB). 



DISCUSSION 

Larval retention and loss were found to depend primarily on larval duration. The 
negative relationship between retention and larval duration is intuitive, in that a longer 
larval duration implies a greater chance of long-distance transport with subsequent loss 
to the system. This is similar to the findings of Leis and Miller (1976), who found that 
larvae of demersal-spawning reef fish (shorter pelagic duration) tended to be found 
closer to shore than larvae from pelagic-spawning reef fishes (longer pelagic duration). 
Some of the proposed physical mechanisms which can transport reef fish larvae back to 
their spawning site operate on the time scale of 2-3 months (e.g., Lobel and Robinson, 
1986), consistent with the relatively high retention found in this study for a 3-month 
larval duration. Late-stage larvae of some reef fishes can occur at great distances from 
suitable adult habitat (e.g., Clarke, 1995; Victor, 1987), but are of unknown importance 
for local population persistence. This issue of long-distance dispersal may, however, be 
important for larval interchange in a metapopulation framework, which will be examined 
elsewhere for insular species in the Flawaiian Archipelago (Kobayashi, in preparation). 
The interannual pattern of retention (Fig. 2 A) is consistent with observed large-scale 
changes in the central Pacific Ocean (e.g., Polovina et al., 1994); however, the exact 
mechanism remains unknown at this point. Larval settlement at Oahu examined by source 
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also appeared to be a proximity effect, with settlement negatively correlated with distance 
from Oahu. Larval duration was an important effect when examining Oahu settlement 
from the adjacent site at Necker (Fig. 6C); however, even at the longest durations, 
the numbers reaching Oahu from Necker did not surpass the number being retained 
around Oahu from Oahu on average (Table 2). Additional modeling is underway to 
better understand connectivity in the Archipelago and will address this on a finer spatial 
resolution, both in terms of source/sink dynamics and oceanographic input data. 

There appears to be an interesting tradeoff between SST and chlorophyll-a with 
respect to seasonal spawning (Figs. 9B, lOB, 12B, 13B). Summer spawning is clearly 
conducive to placing the larvae into higher SST water masses; however, winter spawning 
clearly maximizes chlorophyll-a experienced by larvae. This dilemma does not appear to 
be mediated by seasonal retention (Fig. 2B) or seasonal nonsettlement (Fig. 8B). Early 
summer and late summer may be satisfactory compromises to best optimize these factors, 
thereby keeping both SST and chlorophyll-a at relatively high levels during the pelagic 
stages. The four spawning sites examined in this study generally fall along a latitudinal 
transect, and the resulting site-related patterns in SST and chlorophyll-a are consistent 
with oceanographic work in this area (Polovina et al., 2001 ; Seki et al., 2002). The 
lowest SST and highest chlorophyll-a occur at the northernmost release site of Midway, 
which is well within the TZCF (Transition Zone Chlorophyll Front). At lower latitudes, 
there is a trend for higher SST values, as well as higher chlorophyll-a values. The latter 
may be due to increased nearshore processes (e.g., island effects) enhancing productivity 
around the larger islands in the Archipelago (e.g., Seki et al., 2001). 

In summary, it has been shown that computer simulation may be a useful 
approach towards understanding important aspects of early life history and adult 
spawning ecology. Retention, transport, and environmental variables are shown to be 
expressed in complex spatial and temporal patterns. The utility of this approach depends 
critically on the passivity of larvae. Some late-stage fish larvae have been shown to be 
capable of directional orientation and active movement near the timing of settlement 
(e.g., Kingsford et al., 2002; Leis and Carson-Ewart, 1999; 2000; 2002; 2003; Leis et 
al., 2003; Tolimieri et al, 2004; Jeffs et al, 2003); however, it is quite likely that early 
life history stages (eggs and early-stage larvae) are passive drifters for a large part of 
the pelagic duration, and lobster phyllosoma have very limited swimming abilities. The 
findings of this simulation study can be used to pose further hypotheses and corroborate 
existing empirical evidence. In the latter case, for example, there are observed 
biogeographic patterns in the Hawaiian Archipelago which would benefit from a more 
quantitative mechanistic explanation, such as a higher rate of endemism being found 
at the northerly atolls (DeMartini and Friedlander, 2004), the faunal similarity between 
Johnston Atoll and the Main Hawaiian Islands (Kosaki et al, 1991), and the pattern of 
spread of introduced/invasive species such as the blue-lined snapper Lutjanus kasmira 
(Friedlander et al., 2002). Such corroboration could serve as potential ground-truthing 
for the modeling approach. Additionally, the SST and chlorophyll-a histories provide 
a useful environmental perspective to recent findings emphasizing the importance of 
larval physiological fitness (e.g., Berkeley et al., 2004) towards population maintenance. 
Incorporating demographic variability into the transport-modeling framework is a 
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logical next step. The modeling efforts as described here can help understand and predict 
recruitment success, when coupled with empirical observations and field experiments. A 
better understanding of oceanographic source-sink dynamics and connectivity throughout 
the Archipelago will be helpful towards design of marine protected areas (MPAs) and 
reserves (Cowen, 2002), and will contribute towards more effective management and 
resource utilization in the NWHI. 
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Table 1. Summary of GAM analyses. 



Response variable 


Predictor variable(s) 


Larval retention (all data) 


Year, Month, Site, and 
Duration 


Larval settlement at Oahu (all data) 


Site 


Larval settlement at Oahu (separately by site) 


Year, Month, and Duration 


Larval non-settlement (all data) 


Year, Month, Site, and 
Duration 


Integrated SST history of retained larvae (separately by 
duration) 


Year, Month, and Site 


Integrated chlorophyll-a history of retained larvae 
(separately by duration) 


Year, Month, and Site 



Table 2. Summary of larval settlement at Oahu, aggregated by year and month, from 
different spawning sites and for different larval durations. The average number reaching 
Oahu is out of 5,000 releases; the number of combinations with zero is out of 120 
different year and month combinations per site/duration strata. 



Site 


Duration 


Average no. reaching 
Oahu 


No. combinations with zero 


Midway 


3-month 


0.00 


120 




6-month 


0.00 


120 




12-month 


5.10 


66 


Maro 


3-month 


0.00 


120 




6-month 


14.61 


93 




12-month 


45.15 


14 


Necker 


3-month 


55.61 


82 




6-month 


85.52 


29 




12-month 


98.12 


0 


Oahu 


3-month 


732.73 


8 




6-month 


297.02 


2 




12-month 


126.47 


0 
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shown for illustrative purposes only. 



Addtt.ive compotient. lor retention Addiiivo cornponcinL lor roloniion 

- 0.15 - 0.10 - 0.05 0.0 0.05 0,10 -0 2 - 0.1 0.0 0.1 02 



374 







fifidway Maro Necker Oahu 3 mo. 6 mo. 12 mo. 

Site Duration 

Figure 2. Results of GAM application to larval retention. The predictor variables are year (A), month 
(B), spawning site (C), and larval duration (D). C.I. are ±2 standard errors. 
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Figure 3. Results of GAM application to larval settlement at Oahu. The predictor variable is spawning site. 
C.I. are ±2 standard errors. 
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Figure 4. Results of GAM application to larval settlement at Oahu from Midway spawning. 
The predictor variables are year (A), month (B), and larval duration (C). C.I. are ±2 standard 



errors. 
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Month 




3 mo. 6 mo. 12 mo. 

Duration 



Figure 5. Results of GAM application to larv al settlement at Oahu from Maro 
spawning. The predictor variables are year (A), month (B), and larval duration 
(C). C.L are ±2 standard errors. 
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3 mo. 6 mo. 12 mo. 

Duration 



Figure 6. Results of GAM application to larval settlement at Oahu from 
Necker spawning. The predictor variables are year (A), month (B), and larval 
duration (C). C.I. are ±2 standard errors. 



Additive component for Oahu settlement Additive component for Oahu settlement 



379 




1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 

Year 




3 mo. 6 mo. 12 mo. 

Duration 



Figure 7. Results of GAM application to larval settlement at Oahu from Oahu 
spawning. The predictor variables are year (A), month (B), and larval duration 
(C). C. I. are ±2 standard errors. 



Additive componeitl for noo setOemenl Additive component for non-selliemenl 
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Figure 8. Results of GAM application to larval nonsettlement. The predictor variables are year (A), month 
(B), spawning site (C), and larval duration (D). C.I. are ±2 standard errors. 



Additive component for SST history Additive component for SST history Additive component for SST history 
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Midway Maro Necker Oahu 

Site 



Figure 9. Results of GAM application to SST history of retained larvae after 3-month pelagic duration. The 
predictor variables are year (A), month (B), and spawning site (C). C.I. are ±2 standard errors. 
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Midway Maro Necker Oahu 

Site 



Figure 10. Results of GAM application to SST history of retained larvae after 6-month pelagic duration. 
The predictor variables are year (A), month (B), and spawning site (C). C.l. are ±2 standard errors. 
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Site 



Figure 11. Results of GAM application to SST history of retained larv ae after 12-month pelagic duration. 
The predictor variables are year (A), month (B), and spawning site (C). C.I. are ±2 standard errors. 
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MicJway Maro Necker Oahu 

Site 



Figure 12. Results of GAM application to chlorophyll-a history of retained larvae after 3-month pelagic 
duration. The predictor variables are year (A), month (B), and spawning site (C). C.I. are ±2 standard 
errors. 
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Figure 13. Results of GAM application to chlorophyll-a history of retained larvae after 6-month pelagic 
duration. The predictor variables are year (A), month (B), and spawning site (C). C.I. are ±2 standard 
errors. 
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Figure 14. Results of GAM application to chlorophyll-a history of retained larvae after 12-month pelagic 
duration. The predictor variables are year (A), month (B), and spawning site (C). C.l. are ±2 standard 
errors. 
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